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Abstract. One of the apparent advantages of quantum computers over their classical counterparts 
is their ability to efficiently contract tensor networks. In this article, we study some implications 
of this fact in the case of topological tensor networks. The graph underlying these networks is 
given by the triangulation of a manifold, and the structure of the tensors ensures that the overall 
tensor is independent of the choice of internal triangulation. This leads to quantum algorithms for 
additively approximating certain invariants of triangulated manifolds. We discuss the details of this 
construction in two specific cases. In the first case, we consider triangulated surfaces, where the 
triangle tensor is defined by the multiplication operator of a finite group; the resulting invariant 
has a simple closed-form expression involving the dimensions of the irreducible representations of 
the group and the Euler characteristic of the surface. In the second case, we consider triangulated 
3-manifolds, where the tetrahedral tensor is defined by the so-called Fibonacci anyon model; the 
resulting invariant is the well-known Turaev-Viro invariant of 3-manifolds. 



1. Introduction 

The 1994 discovery of polynomial-time quantum algorithms for factoring and discrete loga- 
rithm [30] motivated a widespread search for more examples where quantum computers outperform 
their classical counterparts. Quantum topology is proving to be an interesting source of results in 
this area. For instance, it is known that quantum computers can simulate certain topological func- 
tors [13] and approximate the value of the Jones Polynomial of a link [2, 1, 14, 28, 31, 35]. In 
this article, we describe another way of applying ideas from quantum topology to quantum algo- 
rithms: approximating invariants of manifolds by contracting tensor networks constructed from 
triangulations. 

Recall that an n-dimensional manifold is a topological space which is locally homeomorphic (i.e., 
topologically equivalent) to M. n . We will consider the cases n = 2 and n = 3: surfaces and 3- 
manifolds. An invariant of n-manifolds (for some fixed n) associates a complex number to every 
n-manifold, such that homeomorphic manifolds are associated to the same number. A simple 
nontrivial manifold invariant is the genus of surfaces, which is computed simply by counting the 
number of handles. If the surface is given in terms of a triangulation, then we can compute the genus 
by the formula (2 - ^vertices + hedges - # faces )/2. The invariants studied in this article are more 
complicated than the genus, but the idea is essentially the same: starting from a triangulation of a 
manifold, we apply some formulas to calculate a complex number; if these formulas satisfy certain 
properties, the number will be an invariant. 

Our goal in the first section is to give an accessible presentation of this idea in the case of 
two-dimensional topological lattice field theories (or TLFTs) attached to finite groups. The idea 
behind two-dimensional TLFTs is based on a theorem of Pachner [25], which states that any two 
triangulations of the same surface differ by a finite sequence of simple combinatorial moves of two 
types. It follows that decorating a triangulated closed surface 1 with tensors which are invariant 
under the two Pachner moves results in a topologically invariant scalar. When the tensors are 



^Unless explicitly stated, surfaces (and more generally, manifolds) are assumed to have no boundary. Requiring 
that a surface be closed ensures that there are also no punctures. 
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defined by the multiplicative structure of a finite group, this scalar has a simple expression in terms 
of the Euler characteristic of the surface and the dimensions of the irreducible representations of 
the group. Mednykh's formula [24] states that this scalar can also be expressed in terms of maps on 
the fundamental group of the surface. The quantum algorithm we present approximates this scalar 
by contracting the TLFT tensor network using well-known methods [5]. Although this algorithm 
does not provide a quantum speedup for the task of calculating surface invariants, it is nonetheless 
useful as a simple prototype of a much more general construction. Moreover, it is conceivable that 
it could provide a speedup for approximating some other quantity involving the relevant group and 
surface. 

In the second section, we study tensor networks arising from triangulations of 3-manifolds. Just 
as surfaces can be triangulated into a collection of triangles glued along their edges, 3-manifolds 
can be triangulated into a collection of tetrahedra glued along their triangular faces. In three- 
dimensional TLFTs, rank-six tensors are assigned to each tetrahedron and attached to one another 
by means of gluing tensors associated to edges of the triangulation. As shown by Turaev and Viro 
[34] , we can choose these tensors so that the three-dimensional versions of the Pachner moves are 
satisfied; this implies that the tensor network of a closed 3-manifold will contract to a topologically 
invariant scalar. This scalar, known as the Turaev- Viro invariant, appears to be difficult to calculate 
classically. No efficient classical approximation algorithms are known, and exact evaluation is NP- 
hard [19, 20]. Approximating it for manifolds specified in terms of a so-called Heegaard splitting is 
a BQP-complete problem [4]. The same problem for so-called mapping tori is DQCl-complete [3]. 
An efficient quantum approximation algorithm is known when the manifold is specified in terms of 
Dehn surgery [16]. In this article, we give an efficient quantum algorithm for approximating the 
Turaev- Viro invariant of a 3-manifold specified by a triangulation. This result is not a consequence 
of the results in [3, 4, 16], as no efficient translation from triangulations to Heegaard splittings, 
mapping tori, or Dehn surgeries is currently known [33]. However, it is possible to efficiently 
translate Heegaard splittings into triangulations, as we will discuss (see, for example, [7].) This 
translation allows us to conclude that approximating the Turaev- Viro invariant of a triangulation 
is BQP-hard, in an appropriate sense. 

As in all of the other aforementioned results, the approximation of our algorithm is additive. 
The approximation scale depends on the quality of the input triangulation as well as the order 
of tensor contraction selected by the algorithm designer. In the worst case, the approximation 
scale is exponential in the size of the input triangulation. On the other hand, we will show that 
every 3-manifold admits triangulations which are efficient in a certain sense, and where a large 
portion of the tensor network can be contracted unitarily, i.e., in a way that does not worsen the 
approximation scale. 

2. Tensor networks and 2D topological lattice field theories 

2.1. Tensors and tensor networks. We begin with a brief review of tensors and tensor networks; 
for a gentler introduction to tensors, see [9]. Throughout this section, V will denote a fixed 
finite-dimensional vector space over C. Recall that an (m,n)-tensor (over V) is an element of 
y®m ^ ^y*^gm^ g uc j 1 a tensor is said to have contravariant rank m, covariant rank n and total 
rank m + n; it can be denoted by 

3132-jn ' 

where each copy of V or V* is associated with a particular index. Contravariant indices are written 
in superscript, while covariant indices are written in subscript. One may think of tensors in a 
basis-dependent way by choosing an orthonormal basis for V. From this point of view, a tensor is 
simply an indexed array of scalars. The above tensor could then be written in Dirac notation as 

Zl M jil 2 'iKl^2,...,im)(il,i2,...,in| , 
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where the sum varies each index over the entire basis of V. We form the product of two tensors by 
composing them as multilinear maps; the use of index notation is crucial here, as two multilinear 
maps can be composed in many ways. For instance, the product of a vector v with a dual vector 
w could either be an inner product v l Wi £ C (if we choose to "act on the same index" ) or an outer 
product v l Wj ZzV ®V* (if we "act on different indices"). Similarly, a matrix A acting on a vector v 
is denoted by Aji^ , the product of A with another matrix B is denoted by AjBj, and the trace of A 
is simply A\. Notice that from the basis-dependent point of view, summation over repeated indices 
is implicit. In the language of tensor analysis, we say that these repeated indices are contracted. 
Clearly, the product of a series of tensors is itself a tensor, and its rank is determined by the indices 
appearing only once in the entire product. In particular, if all the indices are contracted (as in A\), 
then the resulting tensor is simply a scalar in C. 

Tensor networks provide a natural way of depicting and manipulating products of many tensors. 
Imagine an (m, n)-tensor as a widget consisting of one vertex and m + n directed edges, each one 
attached to the vertex at one end. Each edge is labeled by an index of the tensor, and is directed out 
if and only if the corresponding index is contravariant. A product of two tensors then corresponds 
to attaching the two widgets together by appropriately joining the edges corresponding to indices 
shared by the two tensors. A few basic examples of tensor networks are depicted in Figure 1. If 




FIGURE 1. Examples of simple tensor networks (from left to right): a (3,2) tensor 
, the matrix product AjBj, and the matrix inner product A^Bj. 



all of the edges in a tensor network are attached to vertices on both ends, then we say that the 
network is closed. 

Definition 2.1. A closed tensor network (S,M) is a finite directed multigraph 9 equipped with 
tensors M(v), one for each vertex v, so that the edges from v to w in 9 are in one-to-one corre- 
spondence with indices which are both contravariant in M(v) and covariant in M(w). 

A closed tensor network is simply another representation of the product of the tensors M(v), and 
is thus itself a tensor which we will denote by T(9, M). Since each edge is attached to a vertex on 
either end, each index appears twice in T(9, M), which is thus an element of C. This value has a 
simple expression in any given orthonormal basis of V, as follows. A labeling of 9 is an assignment 
of a basis element to every edge of 9- Such a labeling I identifies a particular entry M(y)i G C of 
each tensor M(v). It's then easy to check (e.g., by writing tensors in Dirac notation as above) that 
tensor contraction corresponds to taking sums of products of these scalars. In particular, 

(2.1) T(9,M)= n 

labolings I vertices v 

We now briefly summarize Arad and Landau's quantum algorithm for producing an additive 
approximation of the above number [5]. The algorithm involves a choice of ordering {Mj}™ =1 of 
the vertices of the network, according to which the tensors are contracted one-by-one. At the jth 
stage of the algorithm, we view the network as being divided into three pieces: 
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• the network N(l,j — 1) formed by the tensors Mi, . . . , M,_i; 

• the network N(j + 1, n) formed by the tensors Mj + i, . . . , M n ; 

• the tensor Mj, connected to N(l,j — 1) by some k indices and to N(j + l,n) by some I 
indices. 

The algorithm maintains a register which, at this jth stage, contains an approximation of 2V(1, j— 1), 
viewed ELS £L rank k tensor 2 . We view Mj as a map 

(2.2) : V® k -> T/ 0/ . 

We would like to implement this map on our register, so that it then contains an approximation of 
the rank I tensor N(l,j). To achieve this, we first convert Mj into a square matrix by appropriately 
adding input registers (if k < I) or output registers (if k > I). We then classically compute the 
singular value decomposition of this square matrix. Using a clever trick (see Lemma 3.3 of [5]) we 
can implement the diagonal matrix of singular values as a unitary operator, at an approximation 
scale cost of the operator norm ||M,-|| of Mj, viewed in the form (2.2). At the final stage of the 
algorithm, having contracted all of the tensors, the algorithm approximates the value by means of 
the Hadamard test. The following is a restatement of this result (Theorem 3.4 of [5].) 

Theorem 2.1. Let (3,M) be a closed tensor network where every tensor has degree no more than 
d and is defined over a fixed vector space V . Let {Mj}™ =1 be an ordering of the tensors of (9, M), 
and let A = II?=i ll-^j'll ^ e ^ e product of the norms of the resulting contraction operators. Given 
any e > 0, there exists a quantum algorithm that runs in time poly((dim V) d )n/e 2 and outputs a 
complex number x such that 

Pr(|T(S,M) -x\ > eA) < i 

In the remainder of the paper, we will discuss how this algorithm can be applied to approximate 
certain invariants of triangulated manifolds. 

2.2. Surface invariants from 2D topological lattice field theories. A theorem of Rado from 
1925 states that every orientable surface admits a triangulation [27]. A triangulation is a purely 
combinatorial object, consisting of a list of triangles along with "gluing" rules which tell us how 
to attach the triangles to each other. 3 Each triangle inherits an orientation from the surface, and 
each of its edges are glued to exactly one other triangle. The simplest triangulation of the torus, 
for instance, is formed by making a square sheet out of two identical right-angled triangles, and 
then gluing the opposite edges of the sheet. We point out two ways to change this triangulation 
without changing the topological nature of the torus. First, we could replace the two triangles by 
another pair of triangles, glued along the other diagonal of the square sheet. Second, we could have 
subdivided any one triangle into three by choosing a point somewhere in its interior and adding 
three new edges. These two moves are called the 2-2 move, and the 1-3 move, respectively. It 
is clear that performing any of these moves (or their inverses) does not change the topology of 
the underlying surface, since it amounts to replacing a disk with a disk. A beautiful theorem of 
Pachner [25] states that any two triangulations of the same surface differ by only a finite sequence 
of these moves. Pachner's theorem suggests that we might find topological invariants by looking 
for objects which are unchanged by the 2-2 and 1-3 moves. For instance, it's easy to check that the 

o 

If we ignore the distinction between contravariant and covariant indices, any rank k tensor can be viewed as an 
element of V® k and can thus be stored (up to normalization issues) in a quantum register. 

Note that this is quite different from the notion of "triangulation" used in computer graphics. Our topological 
notion of triangulation does not describe a way of embedding the triangles in ffi 3 , and we allow isotopies for free 
(i.e, triangles are not required to be flat or rigid.) Moreover, we tend to prefer coarse triangulations, since they are 
simpler combinatorial descriptions of the same topological object. 
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Figure 2. Left: a simple triangulation of the torus, and the corresponding tensor 
network. Right: the 2-2 and 1-3 Pachner moves on triangulations. 

quantity ^vertices - #edges + #faces (known as the Euler characteristic) is invariant under these 
moves; in fact, one can show that any linear combination of the number of vertices, edges and faces 
which is invariant under the Pachner moves is a scalar multiple of the Euler characteristic [29] . The 
Euler characteristic of a surface of genus (i.e., number of "handles") g is equal to 2 — 2g. Since the 
genus is a complete invariant of closed orientable surfaces, so is the Euler characteristic. 

Of course, we can imagine defining more complicated algebraic objects which are invariant under 
the Pachner moves; topological lattice field theories (or TLFTs) are based on this idea. A TLFT is 
defined by a complex vector space V and tensors C y ' fc and gij (over V) which satisfy certain tensor 
equations. Given this data, we can associate a tensor to every triangulated surface, as follows. 
Separate out each triangle so that the number of edges is now exactly three times the number of 
triangles. Label each edge of each triangle with an index, i.e., a copy of V . For each triangle, 
assign the tensor C abc , where a, b, and c are the indices of its edges. For each pair of edges s and 
t that should be glued back together to form the original triangulation, assign the tensor g s t- The 
result is a tensor network - the torus case is shown in Figure 2. If an index appears on an interior 
edge of the surface, then it is contracted, as it is covariant in a gluing tensor and contravariant 
in a triangle tensor. The rank of the entire tensor is thus equal to the number of edges on the 
boundary of the initial triangulation. In particular, if the surface has no boundary, the result is a 
scalar. We ensure that this tensor is topologically invariant by requiring that C l i k is invariant under 
cyclic permutations of its indices (as these are the allowed orientation-preserving transformations 
of a triangle), and that both C^ k and satisfy the Pachner moves. For instance, the 2-2 move 
corresponds to the tensor equation 

(2.3) C abc g cd C de S = C ac fg cd C bed . 

This motivates the following definition. 

Definition 2.2. A two-dimensional topological lattice field theory (or 2-D TLFT) is a triple consist- 
ing of a finite- dimensional complex vector space V, a cyclically invariant (3,0)-tensor C % i k (called 
the triangle tensor) and a (0,2) -tensor g^ (called the gluing tensor), both defined over V, which 
satisfy the 2-2 and 1-3 Pachner moves. 

As it turns out, 2-D TLFTs can be classified completely. Notice that the tensor C l i k gki is a 
multiplication map from V (g> V to V, equipping V with an algebra structure. From this point of 
view, equation (2.3) is exactly the associativity property of the algebra product. One can show 
that invariance under the 1-3 move is equivalent to semisimplicity of V. We thus have that two- 
dimensional TLFTs are in one-to-one correspondence with associative, semisimple algebras [15] . A 
natural family of such algebras are the group algebras of finite groups. 

Given a finite group G, the group algebra CG is the vector space of maps from G to C, with 
algebra product given by extending the group product by linearity. The natural action [x ■ f](y) = 
f(x~ 1 y) of G on elements of CG turns the group algebra into a representation of G, called the right 
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regular representation. It's easy to check that the character of this representation, denoted Xreg> 
satisfies 

if a: = 1 
otherwise. 

The tensors of the 2-D TLFT associated to a finite group G are then defined by setting 

rijk _ Xreg(ufc) , _ Xreg(u) 

The tensor C JJlfc is clearly cyclically invariant, and the 2-2 move is satisfied due to the associativity 
of the group product. It is straightforward to check that, with the above scaling, these tensors are 
also invariant under the 1-3 move. 

We can thus associate a number M(G, S) G C to every pair (G, S) where G is a finite group and 
S is a triangulated closed orientable surface, such that M(G, S) = M(G, f(S')) for any choice of 
triangulation S' of S and homeomorphism /. An elegant derivation of this invariant is described 
in [32]; it is given by 

(2.4) M(G, S) = Y, dim(p)^ , 

where x(&) * s the Euler characteristic of S and G is the set of irreducible representations of G. 
Mednykh's Formula [24] states that 

^dim(p)*( 5 ) = |G|^ s )" 1 |Hom( 7 r 1 (S),G)| , 

where if\(S) denotes the fundamental group of S. A beautiful recent proof of this formula, due to 
Snyder, proceeds by computing the value of the TLFT invariant in the group basis (resulting in 
the right-hand side) and in the Fourier basis (resulting in the left-hand side) [32]. 

2.3. Quantum algorithm for approximating the 2D TLFT invariant. The previous section 
described how to create a tensor network associated to a finite group G and a triangulated surface 
S. To apply Theorem 2.1 to this case and produce an algorithm for approximating the invariant 
M(G,S), we need to determine the operator norms of contracting and gij in various ways, 
and discuss the choice of ordering. 

To compute the relevant operator norms, it is useful to write our tensors in Dirac notation. In 
what follows, all sums are taken over the group G. The triangle tensor, viewed as a (3, 0) tensor, 
is equal to 

— v m&>ic> • 



\G\ ^ 

If we write it instead as a (2, 1) tensor and a (1, 2) tensor, then we get 

1 x . 1 



|a)|6)(c| and |a)(6|(c| , 



\ G \ abc=l V\ G \ abc =l 

respectively. It is straightforward to check that the (3, 0) and (0, 3) cases each have operator norm 
v/fCrj while the (2, 1) and (1, 2) cases have operator norm 1. The gluing operator, viewed as a (0, 2) 
tensor and a (1, 1) tensor, is equal to 

£>>|&) and £>)(&|, 

ab=l ab=l 

respectively. It is again straightforward to check that the operator norms are \/\G\ and 1, respec- 
tively. 
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We note that the operators defined above have a rather special structure. Theorem 2.1 is quite 
general and assumes nothing about the operators Mj, implementing them in time polynomial in the 
dimension, which contributes a factor of poly(m d ) to the running time. In our case, for choices of G 
where multiplication can be performed efficiently, the (3, 0) and (0, 3) tensors (scaled appropriately) 
can be implemented in polylog(|G|) time by forming a uniform superposition over G x G and then 
applying the controlled multiplication operator |a)|6)|0) — > |a)|6)|6~ 1 a _1 ). The (2,1) and (1,2) 
tensors can be implemented efficiently in a similar manner. 

The total approximation scale of our algorithm is the product of the operator norms of the 
various tensors, contracted in the order which we choose. This scale will increase by a factor of 
\f\G\ whenever a triangle tensor is contracted as a (0, 3) or a (3, 0) tensor, or when a gluing tensor 
is contracted as a (0, 2) or a (2, 0) tensor. Now, let S be a triangulated orientable surface equipped 
with an ordering {Sj} of the set of all edges and triangles in the triangulation. We will call an edge 
e a cap if the two triangles surrounding e in the triangulation are either both before e, or both 
after e in the ordering {Sj}. Similarly, a triangle t is a cap if the three edges surrounding t in the 
triangulation are either all three before t, or all three after t in the ordering {Sj}. We denote the 
total number of caps by k(S). It is clear that k(S) also denotes the number of \f\G\ operator norm 
factors resulting from contracting the tensor network of S in the order {Sj}. 

By the discussion above, we can now specialize Theorem 2.1 to the case of approximating the 
TLFT invariant M (G, S) attached to an orientable surface S and a finite group G. 

Theorem 2.2. Let G be a finite group and let S be a compact orientable surface equipped with a 
triangulation of size n and an ordering of the set of all edges and triangles; let 

M(G,S) = Y,d X P (S) = IGl^^lHomM^GOl . 
pad 

Then for every e > 0, there exists a quantum algorithm which runs in time polylog(\G\^)n / e 2 and 
outputs a complex number x such that Pr(\M(G, S) — x\ > e\G\ k ^ s ^ 2 ) < 1/4. 

This algorithm is a simple and natural example of applying tensor network contraction to computing 
topological invariants. However, as we now discuss, due to issues of approximation scale and the 
known efficient classical algorithms for computing some of the quantities involved in M(G, S), 
Theorem 2.2 is unlikely to provide any quantum speedup. 

First, it is clear that k(S) is always at least two, since the first and last elements of {Sj} are 
always caps. Moreover, k(S) could be as bad as linear in n; any given triangulation can be worsened 
by inserting the triangulation of the disk shown in Figure 3. If we assume that each edge of the 
disk is contracted in a (1, 1) fashion, then it follows that one of the triangles must be contracted 
(0, 3) or (3, 0). This is simply because the disk has only one exterior edge, and thus only two ways 
(up to direction) of contracting the interior - both of which result in a triangle cap. Inserting such 
a disk will thus increase k(S) by one and the size of the triangulation by three. On the other hand, 
it is easy to come up with triangulations of genus g surfaces with k(S) = 2 (see Figure 3); this 
leads to the best possible approximation scale of \G\. 

There are several quantities in Mednykh's formula that we might try to approximate using 
the above algorithm. First, we could imagine choosing a group for which the dimensions of the 
irreducible representations are known, and then using the above algorithm to approximate x{S)- 
However, since x(S) = 2 — 2g(S) (where g(S) is the genus of S), this approximation scale is 
exponentially large in the genus, even in the ideal case where some p £ G is of dimension nearly 
Moreover, x(S) is trivial to compute exactly using a classical algorithm. 

We might also imagine a scenario (perhaps a black-box problem) in which we know how to 
implement the tensors C l ^ k and gij for some finite group G, but we do not know the dimensions 
of the irreducible representations of G. We could then choose specific inputs of S to try and gain 
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Figure 3. Left: inserting this triangulation of the disk will force a cap, worsening 
the approximation scale by a factor of . Right: triangulations of the torus with 
one and two deleted disks (shaded); gluing g of these together along the boundaries 
in the appropriate way results in a torus of genus g, having only two caps and a 
number of triangles which is linear in g. 



some representation-theoretic information about G. For instance, if g(S) = then M(G,S) = \G\ 
and if g(S) = 1 then M(G, S) = \G\, and so on. Unfortunately, even if we choose the triangulations 
ourselves, the approximation scale is still never better than |G|. 

Finally, we remark that a simple classical probabilistic algorithm for approximating M(G, S) 
has similar performance to the quantum one from Theorem 2.2. In the case of M(G,S), a term 
in the sum (2.1) corresponds to an assignment of group elements to each edge of every triangle 
in the triangulation. An assignment is valid if the cycle of edges around each triangle multiplies 
to the identity and if every glued pair of edges multiply to the identity. Only the valid labelings 
contribute non-zero terms, and each of those terms has the same value: |(^|# tnan gi es -# ed g es . Thus 

M(G,S) = |G| #trianglGS - #edgcs #{valid labelings}. 

This suggests a simple-minded classical algorithm: generate labelings randomly and count the num- 
ber of valid ones among them. A simple analysis shows that this algorithm has similar performance 
to the quantum one. An important distinction is that the quantum algorithm involves a choice of 
ordering of the triangulation. 



3. Three-dimensional topological lattice field theories 

3.1. The Turaev-Viro invariant of 3-manifolds. In 1992, Turaev and Viro introduced a new 
invariant of triangulated 3-manifolds [34]. Their construction, which we now outline, is a three- 
dimensional analogue of the two-dimensional TLFT discussed in the previous section. We begin 
by choosing a finite index set /, a function d : I — > C, and a set J C I 3 of triples from /. These 
pieces of data are referred to by various names in the literature: the elements of / are called labels 
or particle types; the value d{ is called the quantum dimension of particle type i; the elements of 
J are called admissible triples or fusion rules. The final piece of data we need is a function which 
maps each six-tuple (i,j,k,l,m,n) of elements from / to a single complex number: 



i j k 
I m n 



€ C ; 



we will refer to this as the symbol tensor. The symbol tensor corresponds to a tetrahedron whose 
six edges are labeled by (i,j,k,l,m,n) as shown in Figure 4. It will take the value zero unless all 
four of its triangular faces (i, j, k), (k, I, m), (m, n, i), (j, l,n) are admissible; in this case, we will say 
that the six-tuple (i, j, k, I, m, n) is also admissible. The symbol tensor must be invariant under the 




FIGURE 4. Left: a labeled tetrahedron. Middle: the three-dimensional Pachner 
moves. Right: the tensor network associated to two tetrahedra glued together along 
a common triangle; the symbol tensor is denoted by T, and the gluing tensors are 
denoted by g. 



symmetries of the tetrahedron: 
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Using this data, we can turn any triangulated 3-manifold M without boundary into a tensor 
network, as illustrated in Figure 4. Let V, £, and T denote the vertex set, edge set, and tetrahedron 
set of M, respectively. The graph underlying the tensor network will have vertex set £ U T. Two 
vertices in this graph will only be connected if one of them is a tetrahedron t G 7, and the other an 
edge e of t. To turn this graph into a tensor network, we assign the symbol tensor (denoted T in 
the figure) to each element of T, and a gluing tensor (denoted g) to each element of £; this gluing 
tensor will take the value di whenever all of its indices have the same label i, and zero otherwise. 
The gluing tensor g associated to edge e will have rank equal to the number of tetrahedra incident 
on e. 

Every index of the above network will appear in both a gluing tensor and in a symbol tensor, 
so that the overall network will evaluate to a complex number. This number can be specified as 
follows. A labeling of £ is a map <j> : £ — >■ I such that the labels assigned to each triangle are 
admissible. Such a labeling associates to each tetrahedron t£Ta six-tuple (f>(t), consisting of the 
labels of its edges. We denote the value of the symbol tensor of such a labeled tetrahedron by 
\4>(t)\ £ C. As in (2.1), the value of this tensor network is then given by 



e n^)iii^)i- 

labelings <f> e££ t£7 



We would like to use the above to produce a number which depends only on the homeomorphism 
type of the underlying manifold. To achieve this, we must place some additional conditions on 
the symbol tensors and the gluing tensors. These conditions essentially amount to guaranteeing 
that the resulting tensor network is invariant under the three-dimensional Pachner moves, shown 
in Figure 4. First, for any ji, € J such that (ji, j 3 , 34,), (32,j4,h), (ji,33,3e) and (j 2 , 35,3s) 

are admissible, we need 



(3.2) 
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where 5 is the Kronecker delta. Second, for any ji , . . . , jg £ I such that the 6-tuples ( jg , j\ , j'4 , k , j'5 , k ) 
and (j 8 , j 7 , j 9 , j 2 , j 5 , j 3 ) are admissible, we need 
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Third, for any j E /, we need 
(3.4) 
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^ dkdi , 
" 3 k,l: (j,k,l)eJ 

where D = -\J^2j £ j d 2 - is the so-called total quantum dimension. This brings us to Theorem 1.3. A 

of [34]; here, the "initial data" refers to the label set I, the admissible triples J, the dimensions di, 
and the symbol tensor | • |. 

Theorem 3.1. Let M be a triangulated 3-manifold with vertex set V, edge set £, and tetrahedron 
set T. If the initial data specified above satisfies properties (3.2), (3.3), and (3.4), then 

-2|V| 



TV(M) = D~ 



does not depend on the triangulation of M. 



E IKwIIW)! 

labelings <f> eg£ t(L7 



The proof is given in [34]. We remark that Turaev and Viro give a large number of examples of 
initial data which can be used to produce nontrivial invariants by means of Theorem 3.1. These 
examples use the so-called quantum 6j-symbols as the symbol tensors; the initial data, and the 
fact that this data satisfies properties (3.2), (3.3), and (3.4), come from deep facts about the 
representation theory of quantum groups. The simplest of these initial data sets uses only two 
labels (i.e., |/| = 2) and two admissible triples. While this makes for very simple calculations, the 
resulting invariant has a closed form expression involving the Betti numbers of the manifold; these 
can be calculated exactly using standard classical algorithms [18] . 

Although all of our results apply to networks constructed from quantum 6j-symbols as in Turaev 
and Viro's work, as well as more general models considered by Barrett and Westbury [6], we will 
instead write about the so-called Fibonacci model. In this model, there are two labels and three 
admissible triples (up to cyclic permutations). The symbol tensor is very simple to describe, and the 
three essential properties required by Theorem 3.1 can be easily verified. The resulting invariant 
has no known closed form expression, and appears to be computationally difficult for classical 
computers. Indeed, in other settings it is universal for quantum computers [3, 4]. 

3.2. The Fibonacci Turaev- Viro invariant. In the Fibonacci model, the initial data from the 
previous section is quite simple. The label set is / = {0, 1} and the dimensions are 

l + y/l , „ 9 5 + -v/5 



1 , d 1 



and Zr 



2 2 

All triples are admissible so long as they do not have exactly one 1-label. Up to cyclic permutation, 
these are 

J = {(1,1,1), (1,1,0), (0,0,0)} . 

This means that, among all possible six-tuples, there are fifteen admissible ones. Using equalities 
(3.1), one can move between any pair of six-tuples that have the same number of 0-labels. Due to 
this equivalence, we need only define the symbol tensor for the following cases: 



1 1 
1 1 
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Indeed, it is given as follows: 
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recall also that the value must be if (ji,j2, J3, J4,J5, je) is not admissible. Given this initial data, 
we can verify the properties (3.2) and (3.3) directly using a computer program; the property (3.4) 
is quickly verified by hand. We can now specialize Theorem 2.1 to the case of approximating the 
Fibonacci model Turaev-Viro invariant of a triangulated 3-manifold. 

Theorem 3.2. Let M be a triangulated 3-manifold with vertex set V, edge set £ and tetrahedron 
set T. Let d! be the maximum degree (i.e., number of attached tetrahedra) of any edge, and let 
d = m&x{6,d'}. Let {M?}" =1 be an ordering of £ U 7, and let A be the product of the norms of 
the resulting contraction operators. Given any e > 0, there exists a quantum algorithm that runs 
in time poly(2 d )n/e 2 and outputs a complex number x such that 

Pr (J TV(M) -x\> eAD~ 2 ^ < -. 

To analyze the approximation scale of this algorithm, we need to understand the tension between 
A and D~ 2 ^. Recall that during the algorithm of Arad and Landau [5], each time we select a new 
tensor T to contract, we view it as a map 

M T : V® k ->■ V & 

where k is the number of incoming indices, and I is the number of outgoing indices. In our case, V 
is the complex span of / = {0, 1}, i.e., a single qubit. To "contract" T, we convert Mt to a square 
matrix and then implement it approximately as a unitary operator. The cost of this operation is 
a multiplicative increase of A by a factor of ||Mt||, where || • || denotes the operator norm. In our 
case, we can give some simple bounds on ||Mt||. First, if T is a symbol tensor, then 




implies that regardless of k and I, T : |0)® fc i-> |0)®', so that ||T|| > 1. Second, recall that gluing 
tensors take the value if all the indices are equal to i, and zero otherwise. Thus, if T 1 is a gluing 
tensor then regardless of k and I, T : ^ d x so that ||T'|| > \d x \ = (1 + \/5)/2. 

Now suppose that we are building a triangulation by starting from a single tetrahedron and 
attaching tetrahedra one at a time, in the order with which the algorithm will be performed, while 
keeping track of the minimum possible contributions to A' = AD~ 2 ' V L With only the initial 
tetrahedron and its six adjacent gluing tensors, we have 

d\ _ 3 + V5 

Now let us compute A' = Aj^, . Suppose that we are adding the jth tetrahedron by gluing some of 
its triangular faces to the existing triangulation. For our purposes, there are four possible ways to 
attach the new tetrahedron: 

(1) via one of its faces, resulting in three new edges and one new vertex; 

(2) via two of its faces, resulting in one new edge and no new vertices; 

(3) via three of its faces, resulting in no new edges or vertices; 

(4) via four of its faces, resulting in no new edges or vertices. 
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Notice that according to the bounds given above, cases (3) and (4) could possibly result in A'j = 
However, both of these moves reduce the number of triangles in the surface triangulation 
by at least two, while only the (1) move allows us to increase the number of triangles, also by two. 
Since the complete triangulation of M has no boundary triangles, the total number of (3) and (4) 
moves cannot exceed the total number of (1) moves plus one (due to the initial tetrahedron). We 
thus have to assume, for the purposes of best-case analysis of A', that half of the moves are of the 
form (1) or (2). In the case of move (1), we must pay a penalty of at least d\D 2 ps 1.17 while for 
move (2) we pay at least d\ ~ 1.62. In either case, we have IS!- > dfD 2 A'j_ l , and so 

AD -m > 4(^2^/2 = 3W5 m _ 

D 50 

It follows that AZ)~ 2 I V I = 0(2^1). A direct application of Arad and Landau's algorithm to this 
problem thus always results in an approximation scale which is exponential in the number of 
tetrahedra. 

However, a more clever application of the tensor network contraction algorithm could result in 
a much better scale, by taking advantage of the sub-multiplicativity of the operator norm. The 
algorithm designer could carefully choose some adjacent tensors in the network and contract them 
as a single operator, with a possibly lower approximation cost than if he had contracted them one 
by one. 

3.3. Efficient triangulations of 3-manifolds. As we now show, every 3-manifold admits tri- 
angulations such that large portions of the resulting tensor network can be contracted without 
incurring an increase in the approximation scale. Moreover, there are families of triangulated 3- 
manifolds with arbitrarily many tetrahedra, edges, and vertices but constant approximation scale. 
These conclusions come from two observations, which we discuss in detail below: first, that per- 
forming a 2-2 move on a surface by attaching a tetrahedron is a unitary operation, and second, 
that we can efficiently triangulate Heegaard splittings and mapping tori by means of attaching 
tetrahedra to fixed triangulations of handlebodies. 

3.3.1. Attaching a tetrahedron in a 2-2 fashion is unitary. Recall that in Section 3.1, we built a 
tensor network out of a triangulated 3-manifold without boundary by assigning a gluing tensor with 
value di to each edge and the symbol tensor to each tetrahedron, as in Figure 4. We can modify 
this construction slightly to handle 3-manifolds M with nonempty boundary dM. We first assign 
tensors as usual to tetrahedra and edges in M\dM. Then, to edges in dM we assign a gluing tensor 
with value y/dl, with a number of indices equal to the number of tetrahedra adjacent to that edge, 
plus one. The result is a tensor with rank equal to the number of edges in the triangulation of dM. 
We will denote this tensor by T(M). We single out two important properties of this construction. 
The first is that if dM = then by definition we recover the construction presented in Section 3.1 
for 3-manifolds without boundary, and 

T(M) = £> 2 I V I TV(M) . 

The second property is compos ability: if M is a triangulated 3-manifold and N C M a triangulated 
surface such that cutting M along N results in two 3-manifolds M\ and M2, then T(M) is equal 
to the tensor formed by attaching T(Mi) and T(M2) along the indices in N. 

Now suppose M is a triangulated 3-manifold with a boundary dM having k edges. We can think 
of the resulting tensor as a vector |T(M)) € V® k in the complex vector space spanned by labelings 
of the edges of the triangulation of dM by elements of /. Suppose we attach a single tetrahedron 
to M (as in Figure 5), resulting in a new 3-manifold M' . The tensor |T(M')) of M' lives in the 
vector space spanned by labelings of the edges of the triangulation of dM'; this differs from the 
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Figure 5. We can implement 2-2 moves on a triangulated surface (left, grayed) 
by attaching a tetrahedron (middle). Notice that attaching a second tetrahedron 
(right) returns us to the original triangulation. Indeed, by (3.2), the corresponding 
tensor network operator is the identity; attaching one tetrahedron is thus unitary. 



triangulation of dM by a single flipped edge. The map that sends |T(M)) to |T(M')) is known as 
the F-move, and acts trivially on all of the indices except five: 

( 



F : span < 
It is defined by 








f j 




> — > span < 
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32 3l 3 

h h h 




Note that the matrix entries of F^F (i.e., apply the F-move followed by its adjoint, as in Figure 5) 
satisfy property (3.2): 



3 

We conclude that F is unitary. 




J2 Jl J 
33 35 34 
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3.3.2. Efficient triangulations of Heegaard splittings and mapping tori. In the following discussion, 
we will denote the closed, connected, orientable surface of genus g by T, g . Recall that such surfaces 
are completely characterized by their genus. The cylinder S 9 x [0, 1] is a 3-manifold with boundary 
consisting of two copies of E 5 (specifically, the bottom T, g x {0} and the top E s x {1}.) We can 
identify the two ends of this cylinder: 

S g x[0,l] 
(x,0) ~ (x,l) 

to produce a 3-manifold without boundary. Intuitively, the "~" relation means "glue to": each 
point on the bottom of the cylinder is glued to its corresponding point on the top. Now let / be an 
orientation-preserving self-homeomorphism of S g . The mapping cylinder of / is a 3-manifold with 
boundary, defined by 

(S s x [0,1])UE, 



M, 



(x,l)~/(x) 
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where U denotes disjoint union. This construction amounts to "replacing" the top of the cylinder 
S 9 x [0, 1] by gluing on the image of T, g under /; how to glue the points together is specified by 
/. The group of isotopy classes 4 of orientation-preserving self-homeomorphisms of T, g is called 
the mapping class group, and denoted MCG(<?). The simplest examples are MCG(O) = {1} and 
MCG(l) = SL(2,Z). For g > 1, Dehn [12] showed that the MCG is generated by a finite number 
of so-called Dehn twists around simple closed curves on the surface. To perform a Dehn twist, take 
a tubular neighborhood of the curve, cut the surface along the curve, apply a 2ir twist to one end 
of the tube, and then reglue. Lickorish [23] demonstrated a set of canonical curves of size 3g — 1, 
shown in Figure 6. The Dehn twists around these curves generate all of MCG(g). Any mapping 
cylinder M g j can thus be completely specified by giving a word 

/ = /1/2/3 ' ' ' fn 

in these standard generators of MCG(g). We will assume from now on that / is always given in 
this way. 




Figure 6. Left: the Dehn twist around curve c applied to a tubular neighborhood 
of some surface; the image of a curve h under the twist is shown on the right. Right: 
the 3g — 1 Lickorish canonical curves (in bold) on the surface of a handlebody. 

We will consider two ways of turning a mapping cylinder M g j into a closed 3-manifold without 
boundary. The first is the mapping torus, produced by gluing the top of the mapping cylinder to 
the bottom: 

= S g x [0,1] 
9j (s,l)~(/(z),0) ' 

For example, choosing g = 1 and / to be the identity map results in the three-dimensional torus. 
The second construction involves gluing two solid genus-g handlebodies onto M g j, one on top and 
one on bottom. In this case, it might be convenient to visualize the mapping cylinder as a three- 
dimensional annulus, i.e., a shell in the shape of a thickened E g . We can fill the interior of the 
shell with one handlebody, and the exterior with the other handlebody, resulting in a surfaceless 
3-manifold. This manifold is called a Heegaard splitting, and denoted by H g j. Every 3-manifold 
can be specified as a Heegaard splitting for some g and /, but not every 3-manifold can be specified 
as a mapping torus [26] . 

We now show how to triangulate a mapping torus T g j with / = /1/2/3 ■ • ■ fn specified as a word 
in the standard generators. We begin with a triangulation of £ 9 which has a strip of triangles 
along each of the canonical curves shown in Figure 6. This can be accomplished by choosing 
an appropriate triangulation of the genus 1 surface with one deleted disk, and of the genus 2 
surface with two deleted disks, and then gluing g of these together appropriately. To the resulting 
triangulation of S g we attach a shell of tetrahedra which implement f\, that is, a triangulation 

An isotopy between two self-homeomorphisms / and g of a space X is a continuous deformation from / to g 
that maintains the homeomorphism property throughout. Specifically, it is a map H : X X [0, 1] X satisfying 
H(x,0) = f(x) and H(x, 1) = g(x) such that H(-,t) is a homeomorphism from X to itself for every t. 
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of the mapping cylinder M g f 1 . This shell is made up of tetrahedra which are always attached 
to the existing triangulation along two of their triangular faces, resulting in a 2-2 move; how to 
perform a Dehn twist using 2-2 moves is explained below, and illustrated in Figure 7. We repeat 
this procedure for each /j until we have a complete mapping cylinder M g t. To turn this into T g j, 
we need only glue the top triangles to the bottom ones. To turn the resulting triangulation into a 
tensor network, decorate each edge of £ 9 with the gluing tensor di (at an approximation scale cost 
of 0(2 9 )), and decorate each tetrahedron with the unitary 2-2 move tensor from Section 3.3.1 (at 
an approximation scale cost of 1). 

If instead we want to produce a triangulation of H g j from a triangulation of M g j, we need to 
slightly modify the above procedure. First, we start with a handlebody triangulation (instead of 
just a triangulation of S ff ) which also has a strip of triangles along each of the canonical curves on 
its surface. We then attach the triangulation of M g j to the surface of the handlebody. We finish 
by taking another copy of the triangulated handlebody, and gluing the top of M g * to its surface. 
We turn this triangulation into a tensor network just as in the mapping torus case, except we also 
add symbol tensors for each interior tetrahedron of the two handlebodies; the total approximation 
scale is still 0(23). 




Figure 7. We can implement a Dehn twist on a triangulated surface by means 
of a sequence of 2-2 moves. The above annuli represent a strip of triangles living 
on a surface, along a canonical curve. At each step, the dashed edge is flipped. 
Performing the entire sequence above twice performs a full 2n twist on the annulus; 
this is precisely a Dehn twist about the corresponding canonical curve. 

Finally, we specify how to implement a Dehn twist around some loop 7 on a triangulated surface. 
We can implement this twist (in either direction) by performing a particular sequence of 2-2 moves 
on the edges adjacent to 7, as illustrated in Figure 7 (see, for instance, [8] or [21].) We assume that 
each canonical curve will lie along a strip of 2/c-many adjacent triangles for some integer k, as shown. 
By applying an ambient isotopy, we can make the two sides of the strip into concentric circles, with 
the vertices spaced out evenly. Order the vertices {fj}^ =1 along 7 in a counterclockwise manner. 
Now do the following four times: for each i from 1 to A;, flip the clockwisemost edge out of V{ that 
does not lie on 7 (dashed edge in Figure 7). Notice that each flip amounts to shifting the inside 
end of that edge by 2ir/k clockwise along 7, and shifting the outside end by 27r/fc counterclockwise 
along the outer circle of the triangle strip. There are 2k edges, so 4k flips are needed to complete 
a full 2tt twist. Notice that, combinatorially, the resulting triangulation is identical to the previous 
one; one can see this by assigning unique labels to all of the edges before the Dehn twist is applied 
and checking that the incidence matrix of the graph (vertices, edges) is unchanged. This fact allows 
us to compose twists and thus build mapping cylinders M g j for arbitrary Dehn words /. 
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We have thus proved the following. 



Theorem 3.3. There is a constant C such that for every fixed integer g > 0, there exist triangulated 
3-manifolds with an arbitrarily large number of tetrahedra, whose Turaev-Viro invariant can be 
approximated by a quantum algorithm with an additive approximation scale no worse than C 9 . 

Fixing a particular g, we could take the families 

{H gJ : / G MCGG?)} or {T gJ : / G MCG( 5 )} , 

where words are chosen arbitrarily and all manifolds are triangulated according to the above pre- 
scription. Of course, in general we have no guarantee that these families consist of pairwise non- 
homeomorphic manifolds. For instance, H g t = H g j> if / and /' belong to the same double coset 
of the so-called handlebody subgroup of MCG(^), consisting of self-homeomorphisms of T, g which 
extend to the identity on the handlebody of genus g. It's also easy to see that T g j = T gh t h -i 
for any h G MCG(p). However, there does exist an infinite family of pairwise non-homeomorphic 
3-manifolds specified as Heegaard splittings with genus one. These are the so-called lens spaces [26]. 
The fact that such a family can be found in the simplest case g = 1 makes it plausible to conjecture 
that such families also exist for larger g. 

The constant C in Theorem 3.3 is essentially the worst-case norm of contracting the tensor 
network corresponding to a single triangulated "handle", g of which can be attached together to 
form a genus-g handlebody. We can give a bound for C by using the fact that the operator norm of 
a tensor, regardless of how it is contracted, is always bounded above by its Hilbert-Schmidt norm 
(i.e., its L 2 -norm as a vector.) The squared L 2 -norms of the symbol tensor and the gluing tensor 



are 



i,j,k,l,m,n£l 

respectively. We then have 



y i j k =16-4^5 yd. 



5 + Vb 

aj = 

2 



C < (16 - 4V5)*/ 2 (^^\ = (50(45 + 19V5))'/ 2 , 

where t is the minimum number of tetrahedra necessary to triangulate a single handle. While a 
handle can be triangulated with only two tetrahedra, more tetrahedra may be necessary to ensure 
that the canonical curves from Figure 6 all lie along edges of the triangulation. 

3.3.3. BQP-hardness of approximating Turaev-Viro. We now briefly outline a proof that the prob- 
lem of approximating the Turaev-Viro invariant of a triangulated 3-manifold is BQP-hard, in an ap- 
propriate sense. Our proof will be by reduction from the following standard BQP-hard J problem [1]: 
given a quantum circuit U of n gates on g qubits, decide in time poly(n, <?) if |(0® 9 |i7|0® 9 )| 2 < 1/3 
or |(0® s |[/|0® 9 | 2 ) > 2/3. 

Suppose that we are given a quantum circuit U consisting of n gates acting on a total of g 
qubits. We give an efficient procedure for translating such a circuit into a triangulated 3-manifold 
M (equipped with an ordering of the tetrahedra) such that TV(M) approximates |(0® 9 |i/|0® 5 | 2 . 
This translation proceeds in two stages. First, using previous results [4] we can efficiently translate 
U into a word w of length poly(rt, g) in the Dehn twist generators of MCG(g), with the promise 
that TV(H gtW ) is within 1/6 of \{0® 9 |C/|0® 9 | 2 . This procedure relies on an alternative definition of 
TV(M) as the modulus squared of a matrix entry in a certain representation p of MCG(g), and 
the fact that this representation has dense image in the unitary group for all g > 2 [14]. The 
procedure from [4] then guarantees that p(w) approximates U on an appropriate subspace. In the 



^more precisely, PromiseBQP-hard 
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second stage, we use the procedure of the previous section to turn w into an efficient triangulation 
of the 3-manifold H 9jW , along with an ordering of the tetrahedra. It is important that we use a 
small triangulation of the handlebodies, so that the approximation scale factor from Theorem 3.2 
satisfies AD~ 2 '^' < 1. This guarantees that the result of contracting the tensor network of this 
triangulation is still within 1/6 of |(0®^|C/|0®9| 2 . 

The above proof can be viewed as translating quantum states on g qubits into labelings of a 
triangulated surface with genus g and translating local quantum gates into local retriangulations 
of the surface, such that these two translations are consistent. In fact, this idea is used in [21] to 
provide a new quantum error correcting code based on the Turaev-Viro invariant and the associated 
Topological Quantum Field Theory. 

4. Discussion 

A reasonable criticism of the algorithm presented in the previous section is that the approxima- 
tion scale is still exponential in the genus (even though the triangulation could grow unconstrained 
by the genus.) There are a number of reasons why it is unlikely that a significantly better scale can 
be achieved. First, the problem is BQP-hard with the stated scale. Second, in the Heegaard split- 
ting and mapping torus versions of the problem [3, 4], the approximation scales are also exponential 
in the genus. The same is true (with bridge number of the link presentation substituted for the 
genus) for the BQP-complete and DQCl-complete versions of the Jones Polynomial problem [1, 31] 
as well as the Turaev-Viro problem for 3-manifolds specified by Dehn surgery [16]. A recent work of 
Kuperberg shows that for these types of problems, it is unlikely that there are quantum algorithms 
whose approximation scale is not exponentially large in some presentation-dependent quantity [22]. 
We also remark that exact calculation of the Turaev-Viro invariant is NP-hard [19]. 

There are many important similarities between the various results mentioned above. First, both 
the Jones Polynomial and the Turaev-Viro invariant can be viewed as matrix entries (or traces, 
depending on the presentation of the underlying link or manifold) of a representation of a mapping 
class group of a surface. In the case of the Jones Polynomial, the surface in question is the re- 
punctured sphere; in the case of the Turaev-Viro invariant, it is the torus of genus g. From this 
point of view, these invariants are natural candidates for approximation on a quantum computer. 
In a sense, this point of view also makes them natural candidates for BQP-hardness: in both cases, 
the hardness results rely on certain properties of these representations, density in the unitary group 
being the most essential. 

We remark that the BQP results on the Turaev-Viro invariant could be unified by an efficient 
classical algorithm for translating (in any direction) between three presentation types: triangula- 
tion, Heegaard splitting, and Dehn surgery. In the previous section, we described one such efficient 
method which converts Heegaard splittings and mapping tori into triangulations. Similar methods 
may be possible for converting Dehn surgeries into triangulations as well. On the other hand, no 
efficient procedure for converting a triangulation into a Heegaard splitting or a Dehn surgery is 
known at this time [33]. 

It is natural to consider extensions of the above problems in higher dimensions. The Crane- Yetter 
invariant of triangulated 4-manifolds, for instance, is a state-sum invariant analogous to Turaev- 
Viro [11]. The Turaev-Viro invariant is defined by assigning a 6j-symbol to each tetrahedron, with 
one index assigned to each edge; the Crane- Yetter invariant is defined by assigning a 15j-symbol 
to each pentachoron, with one index assigned to each edge and each triangular face. Like the two- 
dimensional TLFT invariant associated to a finite group (but unlike Turaev-Viro) Crane- Yetter is 
known to have a closed-form expression [10] involving the Euler characteristic and the so-called 
signature of the manifold. The Euler characteristic is trivial to calculate exactly from homology; 
the signature has a simple description in terms of the cohomology ring and is thus likely to have 
an efficient exact classical algorithm as well [18]. 

17 



5. Acknowledgements 



G.A. is indebted to Alex Russell and Cris Moore for first telling us about the idea of approxi- 
mating topological invariants by contracting tensor networks. We thank Stephen Jordan, Robert 
Konig and Dylan Thurston for useful conversations. G.A. acknowledges the support of NSERC, 
MITACS, and the U.S. ARO. E.B. thanks Ashwin Nayak and the University of Waterloo Combina- 
torics & Optimization department for accepting him into the summer URA program during which 
his contributions were conducted. E.B. acknowledges the support of NSERC. 

References 

[1] D. Aharonov and I. Arad. The BQP-hardness of approximating the Jones polynomial. ArXiv Quantum Physics 
e-prints, May 2006. 

[2] Dorit Aharonov, Vaughan Jones, and Zeph Landau. A polynomial quantum algorithm for approximating the 
Jones polynomial. In Proceedings of the thirty-eighth annual ACM symposium on Theory of computing, STOC 
'06, pages 427-436, New York, NY, USA, 2006. ACM. 

[3] G. Alagic and S. Jordan. Approximating Turaev-Viro 3-manifold invariants of mapping tori is universal for one 
clean qubit. in preparation. 

[4] Gorjan Alagic, Stephen P. Jordan, Robert Konig, and Ben W. Reichardt. Estimating Turaev-Viro three-manifold 

invariants is universal for quantum computation. Phys. Rev. A, 82(4):040302, Oct 2010. 
[5] Itai Arad and Zeph Landau. Quantum computation and the evaluation of tensor networks. SIAM Journal on 

Computing, 39(7):3089-3121, 2010. 
[6] J. W. Barrett and B. W. Westbury. Invariants of piecewise-linear 3-manifolds. ArXiv High Energy Physics - 

Theory e-prints, November 1993. 
[7] Peter Brinkmann and Saul Schleimer. Computing triangulations of mapping tori of surface homeomorphisms. 

Experimental Mathematics, 10(4):571-581, 2001. 
[8] Richard Brunet, Atsuhiro Nakamoto, and Seiya Negami. Diagonal flips of triangulations on closed surfaces 

preserving specified properties. Journal of Combinatorial Theory, Series B, 68(2) :295 - 309, 1996. 
[9] Michael J. Cloud and Leonid P. Lebedev. Tensor Analysis. World Scientific Publishing, Singapore, 2003. 
[10] L. Crane, L. H. Kauffman, and D. N. Yetter. Evaluating the Crane- Yetter invariant. ArXiv High Energy Physics 

- Theory e-prints, September 1993. 
[11] Louis Crane and David N. Yetter. A categorical construction of 4d topological quantum field theories. In Louis H. 

Kauffman and Randy A. Baadhio, editors, Quantum Topology. World Scientific Publishing, 1993. 
[12] M. Dehn. Die Gruppe der Abbildungsklassen. Acta Mathematica, 69:135-206, 1938. 

[13] Michael H. Freedman, Alexei Kitaev, and Zhenghan Wang. Simulation of topological field theories by quantum 

computers. Communications in Mathematical Physics, 227:587-603. 
[14] Michael H. Freedman, Michael Larsen, and Zhenghan Wang. A modular functor which is universal for quantum 

computation. Communications in Mathematical Physics, 227:605-622, 2002. 
[15] M. Fukuma, S. Hosono, and H. Kawai. Lattice topological field theory in two dimensions. Communications in 

Mathematical Physics, 161:157-175, 1994. 
[16] S. Garnerone, A. Marzuoli, and M. Rasetti. Efficient quantum processing of 3-manifold topological invariants. 

Advances in Theoretical and Mathematical Physics, 13(6): 1601-1652, 2009. 
[17] Allen Hatcher. Algebraic Topology. Cambridge University Press, Cambridge, UK, 2002. 

[18] Volker Kaibel and Marc E. Pfetsch. Some algorithmic problems in polytope theory. In Michael Joswig and Nobuki 

Takayama, editors, Algebra, Geometry, and Software Systems. Springer, 2003. 
[19] R. Kirby and P. Melvin. Evaluations of the 3-manifold invariants of Witten and Reshetikhin-Turaev for sl(2, C). 

In S. K. Donaldson and C. B. Thomas, editors, Geometry of Low-Dimensional Manifolds, Vol. 2, volume 151 of 

Londan Mathematical Society Lecture Notes Series. Cambridge University Press, Cambridge, UK, 1991. 
[20] R. Kirby and P. Melvin. Local surrgery formulas for quantum invariants and the Arf invariant. In C. Gordon 

and Y. Rieck, editors, Proceedings of the Casson Fest, volume 7 of Geomety and Topology Monographs, pages 

213-233. Mathematical Sciences Publishing, 2004. 
[21] Robert Konig, Greg Kuperberg, and Ben W. Reichardt. Quantum computation with Turaev-Viro codes. Annals 

of Physics, 325(12):2707 - 2749, 2010. 
[22] G. Kuperberg. How hard is it to approximate the Jones polynomial? ArXiv e-prints, August 2009. 
[23] W. B. R. Lickorish. A finite set of generators for the homeotopy group of a 2-manifold. In Proceedings of the 

Cambridge Philosophical Society, volume 60 of Proceedings of the Cambridge Philosophical Society, pages 769- 

778, 1964. 

18 



[24] A. D. Mednykh. Determination of the number of nonequivalent coverings over a compact Riemann surface. 

Doklady Akademu Nauk SSSR, 239(2):269-271, 1978. 
[25] Udo Pachner. P.l. homeomorphic manifolds are equivalent by elementary shellings. Eur. J. Comb., 12:129-145, 

February 1991. 

[26] V. V. Prasolov and A. B. Sossinsky. Knots, Links, Braids and 3-Manifolds, volume 154 of Translations of 

Mathematical Monographs. American Mathematical Society, Providence, RI, 1997. 
[27] T. Rado. fiber den Begriff der Riemannschen Flache. Acta Szeged, 2:101-121, 1925. 

[28] M. Rasetti, S. Garnerone, and A. Marzuoli. An efficient quantum algorithm for colored Jones polynomials. 

International Journal of Quantum Information, 6:773-778, 2008. 
[29] J. Roberts. Unusual formulae for the Euler characteristic. ArXiv Mathematics e-prints, January 2002. 
[30] P. W. Shor. Algorithms for quantum computation: discrete logarithms and factoring. In Proceedings of the 35th 

Annual Symposium on Foundations of Computer Science, pages 124-134, Washington, DC, USA, 1994. IEEE 

Computer Society. 

[31] Peter W. Shor and Stephen P. Jordan. Estimating jones polynomials is a complete problem for one clean qubit. 

Quantum Information and Computation, 8:681, February 2008. 
[32] N. Snyder. Mednykh's formula via lattice topological quantum field theories. ArXiv Mathematics e-prints, March 

2007. 

[33] Dylan Thurston, personal communication. 

[34] V. G. Turaev and O. Y. Viro. State sum invariants of 3-manifolds and quantum 6j-symbols. Topology, 31(4):865 
- 902, 1992. 

[35] P. Wocjan and J. Yard. The Jones polynomial: quantum algorithms and applications in quantum complexity 
theory. Quantum Information & Computation, 8:147 - 180, 2008. 

Institute for Quantum Computing, University of Waterloo 
E-mail address: galagic@gmail.com 



19 



